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Abstract. We study Atlas-type models of equity markets with local 
characteristics that depend on both name and rank, and in ways that 
induce a stability of the capital distribution. Ergodic properties and 
rankings of processes are examined with reference to the theory of re- 
flected Brownian motions in polyhedral domains. In the context of such 
models, we discuss properties of various investment strategies, including 
the so-called growth-optimal and universal portfolios. 
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1. Introduction 

In modeling equity market behavior, the goal is to construct models that 
are simple enough to be amenable to mathematical analysis, yet complicated 
enough to capture the salient characteristics of real equity markets. A par- 
ticularly salient characteristic of an equity market is its capital distribution 
curve: 

(1.1) logA; log^(fc)(t) , fc = l,---,n, 

i.e., the logarithms of the individual companies' relative capitalizations (mar- 
ket weights) A*{.)(t) at time t, arranged in descending order, versus the 
logarithms of their respective ranks from the largest company down to the 
smallest. 

The capital distribution curve for the U.S. equity market has shown re- 
markable stability over the last century (see, e.g.. Figure 5.1 of Fernholz 
|12j). and this stability has been captured in the Atlas and first-order mod- 
els introduced in [12j and studied by Banner, Fernholz & Karatzas [3j and 
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others. While Atlas and first-order models are able to reproduce the shape 
and stability of the capital distribution curve, they still fail to provide an 
accurate representation of market behavior. It was shown in [3j that these 
models are ergodic, i.e., each stock spends about the same proportion of 
time in each rank over the long term. While ergodicity may be a nice math- 
ematical property, it does not seem to hold for real markets: in real markets 
the largest stocks seem to retain their status for long periods of time, while 
most stocks never reach the upper echelons of capitalization. Hence, a more 
elaborate model is needed. 

In this paper we generalize the first-order models by introducing name- 
based effects of companies in addition to the rank-based effects that are 
present in the simpler models. The resulting hybrid model has more flexibil- 
ity to describe faithfully the complexity of the whole market; in particular, 
the model has both stability properties and occupation time properties that 
are realistic. 

Relation to Extant Literature: From a different point of view, the Atlas 
model can be seen as a physical particle system, with each company rep- 
resented by a particle diffusing on the positive real line. These individual 
diffusive motions have drift and volatility coefficients that depend on the 
entire configuration of particles at any given moment. Recently Pal & Pit- 
man [20] and Chaterjee & Pal [6], [7] studied such systems, specifically when 
the drift coefficient is a function of the particle's rank and all volatility co- 
efficients are equal to a given constant. Under appropriate conditions on 
the drift coefficients, the system has a unique invariant measure in a lower- 
dimensional space; to wit, the system of the n particles is itself not ergodic, 
but the projected system in a lower-dimensional hyperplane turns out to be 
ergodic, with invariant measure that has an explicit exponential-product- 
form probability density function. Moreover, when the number of particles 
increases to infinity, the system converges weakly to one described by a 
Poisson-Dirichlet distribution on the real line. These considerations are use- 
ful in studying the Atlas model for an equity market, when the volatility 
coefficients are all equal. 

The model is still tractable when its volatility coefficients depend on the 
rankings. Questions of existence and uniqueness for such systems in this 
generality are settled through the theory of martingale problem studied by 
Stroock & Varadhan [22] and Bass & Pardoux [5]. A salient feature is that 
three or more particles may collide with each other at the same time with 
positive probability, or even with probability one, under a suitably "uneven" 
volatility structure. This is a very significant departure from the constant- 
volatility case. Some sufficient conditions on the volatility coefficients for 
the occurrence and for the avoidance of triple (or higher-order) collisions, 
are derived in [16] , by comparison with Bessel processes and with help from 
properties of reflected Brownian motion. 
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The ranked particle system has a deep relation with the theory of multi- 
dimensional reflected Brownian motion studied intensively in the context 
of stochastic network systems by Harrison, Reiman and Williams [13] |14j 
|15j and their collaborators. Recently Dicker & Moriarty [10] provided nec- 
essary and sufficient conditions for the invariant density of semimartingale 
reflected Brownian motions in a two-dimensional wedge to be written as a 
finite sum of terms of product-of-exponential form, by extending the geo- 
metric considerations on the so-called "skew-symmetry" condition. In the 
present paper we use this skew-symmetry condition for the n-dimensional 
reflected Brownian motions to solve the basic adjoint relation introduced in 
a piece-wise constant drift coefficient structure, and find an invariant den- 
sity for the ranked process of the hybrid Atlas model in the form of sum of 
products of exponentials. With this explicit formula we compute the invari- 
ant distribution of the capital distribution curve, as well as the long-term 
average occupation times. 

Another interesting system of ranked particles is Dyson's process of non- 
colliding Brownian motions, which are the ordered eigenvalues of a Brownian 
motion on the space of Hermitian matrices. Recent work by Warren [23] con- 
structs Dyson's process using Doob's /i— transform and Brownian motion in 
the Gelfand-Tsetlin cone, as an extension of Dubedat's work [11] on the rela- 
tion between reflected Brownian motions on the wedge and a Bessel process 
of dimension three. The (infinite) ranked particle systems also appear in 
mean-field spin glass theory of Mathematical Physics. In a further recent 
development, Arguin & Aizenman [Ij analyze robust quasi-stationary com- 
peting particle systems with overlapped hierarchical structures where the 
Poisson-Dirichlet distribution emerges as in [20j. Instead of taking Dyson's 
process or the spin glass theory as our model for rankings in equity markets, 
we obtain the ranked particle system through a general formula of Banner 
& Ghomrasni [3] for continuous semimartingales in the hybrid Atlas model. 

Preview: This paper follows the following structure. We describe our model 
in Section [2l its lower-dimensional ergodic properties in Section [3l the dy- 
namics of its rankings in Section SI its invariant measure and occupation 
times in Section O and some portfolio analysis in its context in Section [6l 
In Section [7] we prove some auxiliary results stated in the main sections. 

Notation: The following notions and notation are useful to describe rankings 
as in [3j. We consider a collection {Q^f^^}i<i^k<n of polyhedral domains in 
M" , where y = {yi, . . . ,yn) G Q^^^ means that the coordinate yi is ranked 
f^th among yi , . . . , y„ , with ties resolved in favor of the lowest index (or 
"name"). Note that for every index i = 1, . . . ,n and rank A: = 1, . . . , n , we 

(i) (i) 

have the partition properties U^^^Q) ^ = M" = D^^^Q)^ ■ 

We shall denote by S„ the symmetric group of permutations of {1, . . . , n} . 
For each permutation p G we consider Tip := n^^iQ^^^^^^ , the poly- 
hedral chamber consisting of all points y G such that yp(k) is ranked 
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k among yi, ...,?/„ , for every k = 1, . . . ,n . The collection of polyhedral 
chambers {7^p}pes„ is a partition of all of M" . 

Since for each y G R" there exists a unique p £ E„ such that y G TZp 
(because of the way ties are resolved), we shall find it useful to define an 
indicator map M" B (xi, • • • , Xn)' = x i— > E S„ such that Xpi^(i) > • • • > 
^p^(n) • III other words, p^{k) is the index of the coordinate in the vector x 
that occupies the A;*^ rank among xi, • • • ,Xn- 

When using matrices and vectors, the vector norm ||x|| := ( Yl^=i ^lY^"^ 
and the inner product (x , y) := Yli=i ^iVi — ^'v ^^^^ x , y G M" , where / 
stands for transposition, are defined in the usual manner. The gradient V 
and the Laplacian A operators on the space of twice continuously differ- 
entiable functions are used in Section El as well as the notations C^(-) and 
C^(-) for the spaces of twice continuously differentiable functions with com- 
pact support, and of twice continuously differentiable bounded functions, 
respectively. 



2. Model 

We shall study an equity market that consists of n assets (stocks) with 
capitalizations X{t) = {Xi{t) , . . . ,X„(t))' which are positive for all times 
< t < oo . The random variable Xi{t) represents the capitalization at 
time t of the asset with index (name) i . 

We shall assume that the log-capitalizations Yi{t) := logXj(t), i = 
1, . . . , n satisfy the system of stochastic differential equations 



(2.1) 



dy^(.t) = i^^gklQ(.)(Y{t))+J, + Jjdt + Y,P^,J dWj{t) 

k=l 3=1 

n 

+ 5^fTaQw(>^(t))dW^i(t), Yi{{))=yi- 0<t<oo 



k=l 



with given initial condition y = {yi , . . . , yn)' ■ The process W{-) := (VFi(-) , 
. . . ,Wn{-))' is an n dimensional Brownian motion. As long as the n di- 
mensional process Y{-) := ( Yi(-) , . . . , Yn{-) )' of log-capitalizations is in the 

polyhedron Q^*^ , the z*'^— coordinate Yi{-) is ranked k^^ among Yi(-) , . . . , Yn{-) 
and behaves like a Brownian motion with drift Qk + li + 1 and variance 
{'^k + 91,1)^ + Pi j ■ The constants 7,7^ and represent a common, 
a name-based, and a rank-based drift (growth rate), respectively; whereas 
the constants ak and pi^j represent rank-based volatilities and name-based 
correlations, respectively. 

Assumption: Throughout this paper we assume that the drift constants 
satisfy the stability condition 

n n 

(2.2) Y,9k + Y.l^ = ^, 

k=l i=l 
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and that the (n x n) matrices 

(2.3) Sp := dmg(ap-i(i) , . . . ,cTp~i(„)) + {pi,j)^<ij<^ 

are positive- definite for every p G S„ , with ak > for every k = 1, . . . , n . 
Equation (j2.ip can be cast in vector form as 

(2.4) dY{t)=G{Y{t))dt + S{Y{t))dW{t), y(0) = y G M" 

for < t < oo , where the functions G : M" ^ M" and 5 : are 

G{y) ■= Yl l^p(2^)" (5'p-i{i)+7i + 7,---,5p-i{n)+7n + 7)' 

P6S„ 

S{y):= l7^p(y)•5p; y G M" . 

Thus (j2.ip is a system of stochastic differential equations with coefficients 
that are piecewise constant, the same in each polyhedral chamber TZp , 
p G . Under the assumption of positive-definiteness in (j2.3p . the sys- 
tem ()2.ip admits a weak solution (Y, W) on a filtered probability space 
(17, .F, {J-'t} ,P) satisfying the usual conditions. By the martingale-problem 
theory of Stroock & Varadhan [23] and the results in Bass & Pardoux [5j, 
this weak solution is unique in the sense of the probability distribution. 



3. Ergodicity 



Thanks to the assumption (j2.2p on the drifts, and taking the average of 
both sides of ()2.ip . we obtain the average log-capitalization process Y{-) := 
Y17=i ^ii') / ''^ in the form 

^ n \ ^ 1 " 

Y{t) = - Vyi+7* + -y"f^fc5fc(i) + - Pi,if^Wi(*), 

1=1 k=l i,i=i 

where Bk{t) := V / l^(.){Y{s))W^{s) ■ k = l,...,n 

for < t < cx) , because of Uf^^Q^*^ = W . Here 5i(-) , . . . , S„(-) are con- 
tinuous local martingales with quadratic (cross-)variations given as {B^ , Bi) 
(t) = 1 6k/ ! and hence are independent standard Brownian motions by 
the F.B. Knight theorem. It follows that the average Y{-) of the log- 
capitalizations Yi{-) , . . . ,Yn{-) grows at a rate determined by the common 
drift 7 , i.e., 

Y(T) 

(3.1) lim — ^ = 7 

holds a.s., by the strong law of large numbers for Brownian motion. Let us 
introduce the column vector i := (1, . . . , 1)' and the subspace 

n := {y G = 0}. 
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Proposition 1. In addition to l\2.S^ - (\2.3\\ . let us impose for every p G E„ 
the following stability condition: 

t 

(3.2) J^(fffe + 7p(fc)) <0; £ = l,...,n-l. 

fc=i 

Then the deviations Y{-) := {Yi{-) - Y {■),... ,Yn{-) -Y{-)) of the log- 
capitalizations Yi(-) , . . . ,Yn{-) from their average are stable in distribu- 
tion: there exists a unique invariant probability measure fi for the H— valued 
Markov process Y{-) , and for any bounded, measurable function / : 11 ^ M 
we have the Strong Law of Large Numbers 

(3.3) lim ^ f f{Y{t))dt = [ /(y)/x(dy), a.s. 
Proof. From (12. ip and (12. 4p . we have 

(3.4) dY{t) = G{Y{t)) dt + S{Y{t))dW{t) ; y(0) = y, 

where y := y - I'y ■ i/n, G{y) := G{y) -7-1 and S{y) := S{y) - 
11' S{y) I n for y £ . By (|2.3p the covariance matrix in (j3.4p is uniformly 
nondegenerate: for all x, y Gil we have 

x'S{y)x = x' S{y)x — x'ii'S'(y) x / n = x S{y) x = '^llp{y) " x'spX 
and 

(3.5) Ao||x|p < x'5(y)x < Ai , 

where Aq (Ai) are the minimum (maximum) of the smallest (largest) eigen- 
values of the positive definite matrices Sp over p G S„ . 

Summation-by-parts, along with (|2.2p and (j3.2p . lead now to 

n n 

y'G{y) = ^yi (s'(p«)-i(i) + 7i) = X] ^P^C^) + 7p«{fc) ) 

i=l k=l 

n n— 1 k 

(3.6) = i/py(„) ^ {gk + ix>y{k)) + X] ^ypHk) - ypyik+i)) X] (^^ + Vw) 

fc=i fe=i e=i 

n 

< cy/n ^ {ypy(k) - ypy{k+i)) < c||y|| < ; y G 11 n 7^p 

A:=l 

where c := n-^/^ maxi<£<„_i ^pg2„ Efc=i(5fc + 7p(fc)) < 0. In the last 
inequality we have used for p G and y G 11 n TZp the properties 
yp(i) > yp(2) > • • • > yp(n) , thus also yp(i) > > yp(„) and 

< n max(y2^i),y2^„)) < n (yp(i) - yp(„))^ 
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Now we consider the one-dimensional process N{t) := f{Y(t)) with 
fin) — (lly|P + ^)^^^ > \\y\\ for yen. An apphcation of Ito's rule gives 

dN{t) = J{Y{t)) dt + [f{y)Y\'S{y)\y^y^^ dW{t) ■ < t < oo , 

I{y) ■■ = {f{y)Y\y'G{y) + ^race(5(y)S(y)')) - 

for y G n . It follows from (j3.5|) . p.6p and the boundedness of S'(-) that 
there exists a constant k > such that /(y) < c / 2 < for > k . The 
diffusion coefficient {f{y)\~'^y'S{y) of A^(-) is a matrix whose entries are 
uniformly bounded by some constants from (13. 5p . 

Thus N[-) is recurrent with respect to the interval (0, k) , and hence 
so is y(-) with respect to i? H 11 for some ball B C M" centered at the 
origin. Combining this with (13. Sp and with the theory of Khas'minskii |18j 
(Theorem 5.1 on page 121), we obtain the existence of a unique invariant 
probability measure for the 11— valued process Y[-) . □ 

The condition (j3.2p ensures that, if yi < 2/2 < • • • < and one subdivides 
at time t = the "cloud" of n particles diffusing on the real line according 
to the dynamics of (j2.ip . into two "subclouds" - one consisting of the ^ 
leftmost, and the other of the n — i rightmost, particles - the two subclouds 
will eventually merge. They will not continue to evolve like separate galaxies, 
that never again make contact with each other; cf. the Remark following 
Theorem 4 in Pal & Pitman [20] for an elaboration of this point in the case 
of the purely rank-based hybrid model with equal variances. 

Corollary 1. Under the assumptions of PropositionUl the long-term aver- 
age occupation time that company i spends in the k^^ rank, i.e., 

1 

(3.7) 0k,i ■= lim — / l (,){X{t))dt, i,k = l,...,n, 

exists almost surely in [0, 1] . 

The resulting array of numbers Ok^i G [0, 1] satisfy X]j=i = YYt=i — 
1 for each "name" i = l,...,n and "rank" k = l,...,n, i.e., 1? : = 
{Gk,i)i<k,i<n a douhly stochastic matrix. Similarly, the average occupa- 
tion time 9p of the market in the polyhedral chamber TZp , namely, 

(3.8) Op := lim i f In.m)) dt 

exists a.s. in [0, 1] for every p € S„ , and we have Ok^i = J2^P' where the 
summation is over the set {p G S„ | p{k) = i} , 1 < i,k < n . 

By Proposition [1] and in particular (|3.3p . the quantity of (|3.7p satisfies 
Ok,= \im ^ Tl „(X(t))dt= hm i ^1 ,.,{Y{t))dt = f,{Q^^), 
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where fi is the unique invariant measure for the process Y(-) of (13.4(1 . Since 

U^^iQf = M" = U^^iQif) , it is obvious that E"=i = E]=i ^kj = 1 
for 1 < i,k < n . Equation (|3.8p is obtained similarly. 



4. Rankings 

Let us now look at the log-capitalizations of the various companies listed 
according to rank, namely 

n 

(4.1) Zfe(t) := k = l,...,n, 0<t<cx). 

i=l 

These are the order statistics Zi{-) > • • • > Zn{-) for the log-capitalizations 
Yi{-), - ■ ■ ,Yn{-) , listed from largest down to smallest. We recall the indicator 
map introduced at the end of Section [H and define the S„— valued index 
process 

^^ ■= p^W = p^W , 0<t<oo, 

so that X(pj(i)(t) > • • • > X(pj(„)(t) . We may thus write Zfc(-) = Y(p,(fc)(-) 
from dH]). 

We shall also introduce the total market capitalization X{-) := Y27=i ^«(') > 
as well as the market weights (relative capitalizations) for the individual 
companies and their ranked counterparts, respectively: 
(4.2) 

Xi(t) (.Zj,{t) 
^J"i{t) ■= -^^^^ , i = l,--- ,n and ^(fc)(t) := ^^^^ , fc = 1, • • • ,n. 

Corollary 2. Under (jg.g|) - (jg.3p and (\3.2^ . the process of ranked devia- 
tions Z{-) := (Zi(-) — yi(-) , . . . , Z„(-) — y(-) )' of the log -capitalizations 
Yi{-),--- ,Yn{-) from their average, is stable in distribution by Proposi- 
tion m and so is the ((M+)"-i X T.n)-valued process (H(-) ,^.) , where 
H(-) := (Zi(-) -Z2(-), • • • , Zn-i{-) - Zn{-))' is the rank-gap process of Y{-) . 

In fact, since Z{-) is obtained by permuting the components of Y{-) ^ 
the stability in distribution of Y{-) implies stability in distribution for Z{-) 
from Proposition [H Moreover, the components of the rank-gap process H(-) 
can be written as linear combinations of those of Z{-) , and the index process 
^. can be seen as ^. = p'^^') , where the range S„ of the mapping p is a 
finite set. Thus, the process (H(-) ) is stable in distribution. 

We shall denote by A^'^{t) := Kz^~Zj{t) the local time accumulated at 
the origin by the nonnegative semimartingale Z]^{-) — Zj{-) up to time t for 
1 < fc < j < n, and set A°'i(-) = = A"''^+i(.) . Then from Theorem 2.5 
of Banner &: Ghomrasni [4j it can be shown that we have for fc = 1, . . . , n , 
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< t < oo the dynamics 

n 

dZfc(t) = ^l « {Y{t))dY,{t) 

(4-3) 

e=k+i i=\ 

Here (t) is the cardinahty of the set of indices of those random variables 
among Yi{t), • • • , Yn{t) which have the same value as Zk{t) , i.e., Nk{t) := 
\{i : Yi{t) = Zj,{t) }| . Note that under the assumptions on the coefficients, it 
can be shown that the finite variation part of the continuous semimartingale 
Y{-) in (|2.1|) is absolutely continuous with respect to Lebesgue measure a.s., 
and it follows from an application of Fubini's theorem and an estimate of 
Krylov [E] that the Lebesgue measure of the set {t : Yi{t) = Yj{t)} is zero 
a.s. for 1 < i, j < n . Thus we can verify the sufficient conditions (2.11-12) 
of Theorem 2.5 in [4]. 

Each local time A'^'^(-) is flat away from the set {0 < t < oo | Zk{t) = 
■ ■ ■ = Zi^t)} ; it grows only when the corresponding coordinate processes 
collide with each other. Examples in [5], [T^ study such multiple collisions 
of order three or higher, and use comparisons with one-dimensional Bessel 
processes in a crucial manner. Here again, the nonnegative semimartingale 
Zk{-) — Zi{-) is compared to an appropriate Bessel process. Since a Bessel 
process with dimension 5 > 1 does not accumulate any local time at the 
origin (a consequence of Proposition XI. 1.5 of [21j), appropriate comparison 
arguments yield the following result; its proof is in Appendix l7.1[ 

Lemma 1. Under the local times A^'^(-) generated by triple or higher- 

order collisions are identically equal to zero, i.e., A'^'^(-) = for 1 < k,£ < 
n and \k — l\ > 2 , and ( |^.3D takes for k = 1, . . . ,n , < t < oo the form 

(4.4) dZfc(i) = E ^Qf^^^ dY,{t) + -(^dA'^^'^+\t) - dA^-'^^t)) . 

i=l 

Proposition 2. Under the assumptions (\2.2^ . and (jg.gl) . we obtain 

the following long-term growth relations, in addition to those of l\3.1^ : all 
log-capitalizations grow at the same rate 

,4.5) ,.1,...,„ 

almost surely. This holds also for the total market capitalization 

1 1 " 

(4.6) hm -logX(r) = hm -log ( Vx,(r)) =7, a.s. , 

i=l 

thus the model is coherent: in the notation of ([73 have 



1 

(4.7) lim —log ;Ui(T) = 0, a.s.; z = l,...,n. 
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Proof. It follows from Corollary [2] that 

lim ^{Zk{T) - Zk+i{T)) = , a.s. ; k = l,...,n-l. 
Combining this with (|2.ip . (j3.7p and (j4.4p . we observe 
lini ^(Afc-i.^'(T) + A'=+i''=+2(r) - 2A^'^'+i(r)) 

T — *CXD 21 

n n 

= 5fc + X] ^^^'^'^ ~ {pk+i + ^ 7j^fc+i,i) = Qk- Qk+i , a.s. 
j=i 1=1 

where := gk+Y17=i li^k,i for A; = 1, ... , n— 1 . Adding up these equations 
over k = i , . . . ,n — 1 yields 

(4.8) hm ^(A^-i'^(r) - A^'^+i(r) - A"-i'"(r)) = 5, - 0n , a.s. 



for each i = 1, . . . , n ; and adding up over all these values of ^ , we obtain 
(4.9) lim -LA"-i,n(r) = 

In conjunction with (14. 8|) . we obtain from (14.90 that for A: = 1, . . . , n : 

1 " 

hm — {K^-^'\T) - A'=''=+i(r)) = 0,. = + V 7,^,,, , a.s. 

T — s-oo 11 ^ — ' 

i=l 

From this, (j4.4p . and the strong law of large numbers for Brownian motion, 
we get the long-term average growth rate of ranked log-capitalizations: 

hm — — — =7, a.s.; K = l,...,n. 

This yields (I4.5|l . the elementary inequality exp{yp!/(i)} < X^^Li expjyj} < 
n exp{2/p!/(i)} for y £ M" implies (|4.6p . and equation (j4.7p is a direct con- 
sequence of (|4.5p and (|4.6p . □ 

Corollary 3. Under l\2.l^ . (1^.51) and (jg.^p . t/ie long-term average occupation 
times Ok,i of <\3. ?P satisfy the equilibrium identity 



(4.10) X] ^'^'^^'^ + = ° ' « = l>--->^- 

fc=i 

Indeed, by substituting (j4.5p into (j2.ip we obtain the a.s. identities 

n „j' 

(i){Y{t)) dt = —ji a.s.; i = l,...,n 



1 /" 
lim -^gk 1 (.) 



and so in conjunction with (|3.7p we deduce (j4.10p . 

Example 1. Suppose that the rank-based growth parameters are given as 
gn = (n- — l)g , gi = ■ ■ ■ = gn-i = —g < for some g > Q . This is the 
"Atlas configuration", in which the company at the lowest capitalization 
rank provides all the growth (or support, as with the Titan of mythical lore) 
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for the entire structure. Suppose also that the name-based growth rates 
71, • • • ,7n satisfy YJi=i 7* = and maxi<j<„7j < g . 

It is then easily checked that (|2.2p and (13. 2p are satisfied. By Corollary [U 
the average occupation times {9k,i\ exist a.s. We shall provide an explicit 
expression for the O^^i under an additional condition (15. 7p on the correlation 
structure, in Section [5.21 For the time being, let us just remark that in this 
case we get directly from (I4.10p the long-term proportions of time 

On,i = -fl - — ) ; i = 1,... ,n 
n V g J 

with which the various companies occupy the lowest ("Atlas") rank. □ 



5. Invariant Measure 



5.1. Reflected Brownian Motions. Observe now from ()4.4p the repre- 
sentation for the vector S(-) = (Hi(-), . . . of gaps in the ranked 
log-capitalizations Hfc(-) := Zk{-) — , k = 1, - ■ ■ ,n— 1: 



(5.1) m = 

Here we have set C,{-) :- 



■ (0) + C(t) + 9^A(t) ; < t < oo 
(Ci(-),...,Cn-i(-))' with 



l^,^{Y{s))dY{s) 



/ V) {Y{s))dY{s)- 



and we have introduced the vector A(-) := (A^'2(-), . . . , A"'~^'"(-))' = (Asi(-), 
. . . , A=,^_j(-))' of local times, as well as the ((n — 1) x (n — 1)) matrix 



(5.2) 



( 



\ 



1 

-1/2 



-1/2 
1 



-1/2 
-1/2 



1 

-1/2 



-1/2 
1 



This rank-gap process H(-) in (|5.ip belongs to a class of processes which 
Harrison & Williams [11] , [15] and Williams [25] call "reflected (or regulated) 
Brownian motions" (RBM) in polyhedral domains. 

The process H(-) has state-space (]R_|_)"~"^ and behaves like the (n — 
1)— dimensional continuous semimartingale (^(•) on the interior of (M^)""^ . 
When the face := {(zi, . . . , z„_i)' G (]R+)"-i | Zfc = 0}, A; = l,...,n-1 
of the boundary is hit, the A;*'' component of A(-) increases, which causes 
an instantaneous displacement (reflection) in a continuous fashion. The 
directions of this reflection are given by the entries in , the k^^ column 
of the matrix 9^. For every principal submatrix 9^ of 9^, there exists 
a non-zero vector y such that 91 y > , and so the reflection matrix 91 
satisfies the so-called completely— S (or "strictly semi- monotone" , see Dai 
& Williams [9] for details) condition for 5 = ( 



\n-l 



12 T. ICHIBA, V. PAPATHANAKOS , A. BANNER, I. KARATZASs AND R. FERNHOLZ 

Let US define the differential operators A and Vk , acting on C^((M+)"^^) 
functions: 

(5.3) 2,% dzkdz, ^ dzk 

[Vkf]{z):={xk,Vf{z)); k = l,...,n-l, z G (M+)"-i . 
Here (aA;/(-))i<fc/<n-i is the corresponding covariance matrix with entries 
afc/(p) := {(rl + ■ lik=£} - <^k ' l{fc=f+i} - c^fc+i • l{fc=<?-i} 

n 

/C.N + E (Pp(k),m - Pp{k+l),m)iPp{e),m " Pp(f+l),m) 

^ ' ' m=l 

+ E {'^<^(Pp{l3),a - Pp{f3+l),a) + <7a+l (Pp(/3+l),a+l ~ Pp{/3),Q+l)} 

for k,£ = l,...,n — l,p E S„ ; whereas the ((n — 1) x 1) vector is the k^^ 
column of the reflection matrix 9^ . We also define the ((n — 1) x 1) drift 
coefficient vector 6(-) := (&i(-), . . . , &,„_i(-))' for the semimartingale C(') > 
with components 

(5.5) &fc(p) := c/fe+7p-i(fc)-5fc+i-7p-i(fc+i) ; A; = 1, . . . , n-1 , p S S„ . 

From Corollary [5] we know that there exists an invariant measure z/(- , •) 
for the ((M+)"~-'^ x valued process (H(-),*p.) . Let us denote by fo(-) 
the marginal invariant distribution of H(-) . As a consequence of Ito's for- 
mula and the formulation of the submartingale problem studied by Stroock 
& Varadhan [22j and Harrison & Williams [14j, we obtain a characterization 
of the invariant distribution z^(-,-) for (H(-),^p.). 

Lemma 2. Assume !\2.2^ . (jg.gp and For each k = 1, . . . ,n-l there is 

a finite measure z^ofc(") > absolutely continuous with respect to Lebesgue mea- 
sure on the k^^ face > such that the so-called Basic Adjoint Relationship 
(BAR) holds for any C^— function f : (M+)"~^ M, namely 

(5.6) / [^/](z,p)dz.(z,p) + iV / [Dkf]iz)d,.ok{z)=0. 

This condition is necessary for the stationarity of z^(-,-). A proof of 
Lemma [2] is given in Appendix 17.21 It is not easy to solve (j5.6p in gen- 
eral; however, following Harrison & Williams [15], we may obtain an ex- 
plicit formula for the invariant joint distribution z^(-, •) under the so-called 
sicew symmetry condition between the covariance and reflection matrices; 
see Proposition [3] and Corollaries H] & [5l 

Lemma 3. Assume that the rank-based variances {o"^} grow linearly, and 
that there are no name-based correlations in l\2.3\i . i.e., 

(5.7) al- al = al- al = ■ ■ ■ = al_i - al , pij = ; l<i,j <n. 
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Then the components of covariance matrix 21 = (ci/c,^)i<,t £<n-i /'"^"^ ( |5.^P 
become 

o"fc + • l{fc=^} - • l{fc=£+i} - o-fe+i • ^{k=e-i} 

and does not depend on the permutation p G Sn • Moreover, it satisfies the 
so-called skew symmetry condition 

(5.8) (22)-fiD-2:)^)-22t)^^ = 0; l<fc,^<ri-l. 

Here we have introduced the diagonal matrix D := diag{'QV) , and the ((n — 
1) X (n — 1)) matrix ^ := / — 91 from the reflection matrix 9^ in (j5.^ . 

Lemma [3] is proved by straightforward computation; details are in section 
5.5 of [T6]. Note that under (|5.7p the operator (|5.3p still depends on the 
permutation p through the drift component 6(p) for p G S„ in (j5.5p . 



Proposition 3. Under (fOjl . (fOj) . (fXlj) and ( [577D , t/ie invariant joint 
distribution v{-) of the x Tin)— valued process (^(O;^-) ^-^ 

(5.9) i/(A X i?) := [Y^Il ^qi ) E / ( - ' ^)) ' 

/or any p € S„ where Ap := (Ap^i, . . . , Ap^„._i)' is the vector with compo- 
nents 

(5.10) Ap,, := ^ii%^i±^ ; p E E„ , 1 < ^ < n - 1 , 

for any measurable sets A C (M+)"~"'^ and B C S„ . /n particular, the 
density p(-) o/ the marginal invariant distribution voi') of H(-) has the 
sum-of-products-of-exponenentials form 

n— 1 

(5.11) P(^):=( E n Vfc) E «^P(-^V,^)); ^e(M+)"-^ 

Proof. First, we carry out a linear transformation of the state space to re- 
move the correlation between the components of H(-) ; this is possible, since 
the covariance matrix 21 does not depend on the index process *p. , un- 
der (j5.7p from Lemma [3l Let il be the unitary matrix whose columns are 
the orthogonal eigenvectors of the covariance 21 , and let £ be the cor- 
responding diagonal matrix of eigenvalues such that £ = il'2lil. Define 
S(-) := £~^/^itH(-) . By this deterministic rotation and scaling, we obtain 

(5.12) H(t) = H(0) +C(t) +^A(t) ; < t < cxd 

from where C(-) = £~^/^HC(-) is a Brownian motion with drift coeffi- 
cient b{-) := £~^/^il6(-) where b{-) is defined in ()5.5p . We may regard H(-) 
as a refiected Brownian motion in a new state space © := £~^/^il (M+)"~"'^ 
with faces '■= £^^^^H5'fc , /c = l,...,n — 1. The transformed reflection 
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matrix := H^^/^^IH can be written m = {m + £l)C = (ti, . . . ,^n-i) , 
where e := ^-^/^ , S) = diag(2t) , 91 = = (ni, . . . , Hn-i) , := 

— ^ = (qi, . . . ,^n-i) • The constant vectors rfe,'qfc,nfc , k = 
1, . . . , n — 1 are {{n — 1) x 1) column vectors. 

The corresponding differential operators A , "Dk and their adjoints A* , 
are defined by 

[Af] {z, p) := J A/(z) + (Kp) , V/(z)) , [Vkf] (z) := , V/(z)) , 
(5-13) ^ ? ^ 

[A* f] {z, p) := - A/(z) - (6(p) , V/(z)) , [Dlf] (z) := (F, , V/(z)) , 

where := nfe-qfc + (nfc, qfc)nfc for A; = 1, . . . ,n-l , z G (M+)""^ , p G S„ . 

With these differential operators as in Lemma [21 we obtain the (BAR) for 
the process (H(-),*p. ) and its invariant distribution ?(•,•); i.e., for every 
fe = 1, . . . ,n — 1 , there exist a finite measure {^'ofc(")} which is absolutely 
continuous with respect to the (n — 2)— dimensional Lebesgue measure on 
aiid such that for any C^— function / : (3 ^ M we have 

n— 1 

(5.14) / [Af]iz,p)dViz,p) + ly2 f [Vkf]iz)dVokiz)=0. 

Our argument, especially from here onward, relies heavily on the elaborate 
analysis given by Harrison &; Williams [13], [IS]- The main distinction be- 
tween their setting and ours, is in the drift coefficient 6(-) , which here varies 
from chamber to chamber as well as within each chamber, and is evaluated 
along the path of the index process *p. . Here, however, we can use the 
following. 

Lemma 4. The following two conditions are equivalent : 

(i) For each collection of constants {gk,^i] 1 ^ i^k < n} , there are 
(n — 1) — dimensional vectors Ap := (Ap_i, . . . , Ap^„_i)' for p G S„ , 
such that a probability measure in the form of sum of products of 
exponentials 




for measurable sets ^ C S and C S„ , satisfies ( |5. j^p for /(•) G 
C^(S) , where c in (\5.15\i is a normalizing constant. 
(ii) The covariance and the direction of reflection satisfy the skew sym- 
metry condition l\5.8\i . 
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Indeed, substituting ()5.15|) into (I5.14p and combining the summation over 
p e $]„ , we observe that 



holds for / G C|(S) , where the expression in the curly bracket corresponds 
exactly to the BAR condition studied in [15j with some differences in nota- 
tion. This way, we may reduce our problem to the case of [15]. Following 
the proof of Lemma 7.1 in [15J, we observe that the condition (i) in Lemma 
m above is equivalent to the following conditions (iii)-(iv), where: 

(iii) [A*p.] (•, •) = in e X E„ , and 

(iv) [P*p.] (•, •) = 26fc(-)p.(-,J on §fc X for A; = 1, . . . , n - 1 . 

Here the adjoint operators A* , are defined in ()5.13p . 

Then the same reasoning as in the proof of Theorem 2.1 in [15] yields 
our Lemma m and we obtain Ap = 2(7 — OT£3)~^6(p) for p G S„ along 
the way. This gives the invariant distribution of H(-) in (I5.12|) . Now 
transforming back to H(-) , we obtain (IS.lOp . (15. 9p and then ()5.1ip . □ 

Example 2. With 7^ = , pij = , 1 < i, j < n and = • • • = o"^ , we 
recover the case studied by Banner, Fernholz &; Karatzas [3] and Pitman &; 
Pal [20]. Our Proposition [3] is an extension of their results. □ 

5.2. Average Occupation Times. The long-term average occupation time 
9p of the vector process X(-) in the polyhedral chamber TZp of (|3.8p is the 
probability mass z^i(p) := i^((M+)"~"'^, p) assigned to such a particular 
chamber by the marginal invariant distribution of the index process , 
which we can compute directly from (j5.9p for p G . 

Corollary 4. Under the assumptions of Proposition\^ the long-term aver- 
age occupation time 9p of X(-) in the chamber TZp for p G Sn ; and the 
long-term proportion 9k^i of time spent by company i in the k^^ rank as in 
^3.7\) , are explicitly given by the respective formulae 




/ [Af]{z,p).pp{z)dz + -Y I [Vkf]{z)Ypp{z)dz 






n-1 



n-1 



(5.16) 



and 




1 



Here Ap is in (jJ.iOp . and the summation for Oj^i is taken over the set 
{p G S„ I p(A;) = i } for 1 < i,k < n . 



From Corollary [31 the average occupation times {Ok,i) satisfy the equilib- 
rium identity (j4.10p . As a sanity check, we verify this for the special case 
(|5.16p through some algebraic computations in Appendix 17.31 
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Example 3. It should be noted that in the presence of name-based vari- 
ances, (j5.16p can fail significantly. Consider the case where n = 3, with 
7i = 0, for i = 1, 2, 3; fjfc = o" > 0, for A; = 1, 2, 3; (73 = (7 > 0, (72 = 0, and 
gi = —g; all pij is zero for i,j = 1,2,3 except ^3^3 = a. In this case, 
Yi(-) and l2(") will vibrate quietly in the middle with variance rate a'^, 
while Y^{-) , with much greater variance rate (cr -|- p)^ , will be wandering 
far and wide. From Corollary [1] and ()4.10p we obtain 



(5.17) 1} = {ek,^) 



l<i,k<3 



1-a 


1-a 


2 


2 


a 


a 


1-a 


1-a 


2 


2 




where the parameter a is in the interval (1/3,1/2) for p > . The upper 
bound 1/2 is obtained as limp_»oo ^1,3 • Without name-based variances, 
i.e., if the pij were all zero, the l^(-) would each spend the same proportion 
of time in every rank. This gives the lower bound 1 / 3 . □ 

Example 4. Let us consider a numerical computation of {9k i) for descend- 
ing name-based drifts 7^ and ascending rank-based drifts gk , e.g. n = 10 
and al = 1 + k, as well as gk = —1 for A: = 1,... ,9, gio = 9, 7i = 
1 — ( 2i ) / (n -|- 1) for i = 1, . . . , n . This is a rather extreme case of Exam- 
ple [H with g = 1. The overall maximum is 61^1 = 0.5184, and the overall 
minimum is ^1,10 = 0.00485 . The company "i = 1 " stays at the first rank 
longer than any other companies, because of its relatively strong name-based 
drift; whereas the company "i = 10" stays at the first rank only for a tiny 
amount of time, because of its relatively poor name-based drift. 

Figure [1] shows a gray scale heat map for the different values of j of 

course we know from Example [1] that 0io,i = i /55 , i = 1, • • • , 10 . □ 

For a larger number of companies, say n ~ 5000 , it seems rather hopeless 
for the current computational environment to perform direct computations 
of Ok^i via the sum of (|5.16p over (n — 1)! permutations in general. 

5.3. Capital Distribution Curve. The capital distribution curve is the 
log- log plot of market weights in descending order, as in (jl.ip . The em- 
pirical capital distribution curves, for the U.S. stock market 1929-1999, are 
shown in ^12\ (Figure 5.1 on page 95). Our next result computes the capital 
distribution curves directly from Proposition [3l using change of variables. 

Corollary 5. Under the assumptions of Proposition 0, the ranked market 
weights //(i)(-),... ,Pi^n)i') (|i.-?p . I \4-^ have invariant distribution with 



n—l n 



-1 



(5.18) p(mi, . . . ,m„_i) = ^ [^p " 11 '^P''^ " ( 11 

p6S„ k=l j=l 

as its density, for < m„ < nin-i < • • • < rni < 1 and nin = 1 — mi — 
• ■ ■ — mn-i ■ Here we set Ap^o = = Ap^^ , P & for notational simplicity. 



HYBRID ATLAS MODELS 17 

Moreover, the log-ranked market weights Cfc(-) := log/x(;,)(-) have invari- 
ant distribution with density 

n-l 

(5.19) p(ci, . . . ,c„_i) = [^p • n (Vj- • e-^^P'^-^p^^+i)'^^) • e^''--!^" 

p6S„ j=l 

/or -oo < c„ <•••< C2 < ci < , c„ = log (l - J2]Zl e^') ■ 

From the invariant density functions given by (15.111) . ()5.18p and ()5.19p . 
the piecewise hnear capital distribution curve p.l|) has the expected slope 

2m F'^ r iog/^(fc+i) -logm i ^ ^''(-fc) ^ Epgs„ ^pAp,fc 

^' ^ L log(A: + l)-log/t J log(l + A:-i) log(l + A;-i) 

between the k^^ and the {k + 1)*^* ranked stocks for k = 1, . . . , n — 1 , and 
the initial value 



E^(log/i(i))=E'^(ci)=E^ 



log l+e"=i+e"^=i+'=2;_^. . ._^g 



(Hi+=.2), . . ._i_p-{=iH l-H,i_i) 



for the first rank. From ()5.9p this expected initial value may be obtained 
through a Monte Carlo simulation of generating (n — 1) independent expo- 
nential random variables with intensities Apj- for j = l,...,n— 1, pGS„. 
From (I5.20p we obtain the following simple criterion for convexity (or con- 
cavity) of the expected capital distribution curves. 

Corollary 6. Under the assumptions of Proposition\^ a sufficient condition 
for the expected capital distribution curve log A; i— > E'^(log /i(fc)) under the 
invariant distribution v to be convex (respectively, concave), is that 

(5.21) Ap,fc+ilog(l + ^) -Ap,fclog(l + ^) >0; V p E S„ 

(resp., < ) hold on each interval [log k, \og{k -\- 2)] for each k = 1, ... n — 2 , 
where Xp^k is given in (j5. jOp . 

Example 5. Let us consider the first-order Atlas model which is a combina- 
tion of the "Atlas configuration" in Examples [1] with the further restrictions 
of Example [21 to wit, gn = {n — l)g , gi = • • • = gn~i = —g < for some 
5 > , as well as 7^ = , pij = , I < i,j < n and erf = • • • = o"^ = 
> for some o"^ > . From Corollary [6l the expected capital distri- 
bution curve is convex but almost linear for larger k . Indeed, the quantity 
Ap,felog(l -|- A;^^) = 2(gk / cr^) • log(l -|- A;^^) increases in A; > 1 , and con- 
verges to one, as A | 00 , for all p G S„ , and so the difference in (j5.2ip is 
positive for each k = 1, . . . , n — 2 but decreases to zero quite rapidly in the 
order of 0{k~'^) , as A | 00 . Another explanation of such linearity ("Pareto 
line" ) of the capital distribution curves from an application of Poisson point 
processes can be found in Example 5.1.1 on page 94 of [12] . □ 

Example 6. Suppose now that we change only the rank-based variances in 
Example El namely, we take linearly growing variances af. = ka"^ for some 
^2 > , A = 1, . . . , n . Then Ap,fc log(l A"^) = (4Ac/ / [(2A + 1)0"^]) ■ log(l + 



dYiit) 
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k~^) is decreasing in > 1 for every p G S„ , and so the difference in 
(j5.2ip is negative for each k = l,...,n — 2. Thus, from Corollary [6], the 
expected capital distribution curve becomes concave. □ 

Example 7. Consider a "pure" hybrid market defined by 

^--fidt + pi dWi (t) if Yi{t)^ y(„) (t) , 

ig-ji)dt + pi dWiit) if Yi{t) = Y(n)it), 

for i = 1, . . . , n and t >0, where > 0, pi > 0, and g = J2^=i li- We claim 
that for large n and g ^ the capital distribution curve for this market is 
convex. 

For large n and g » 7j, the process [Yi{t) — y(„))(t) is approximately 
exponentially distributed for t outside the set where Yi{t) = y(„)(t), with 

p{yi(^)-y(„)(^)>x}-e-°'^ 

where Ui = pf/2^i. It would seem that for large enough n, the Atlas stock 
would perform a role similar to a local time process, reflecting the stocks 
away from the Atlas position. 

Now, if we let Y represent a generic member of the market, the distribu- 
tion of Y(t) — Y(„)(t) will be mixed exponential, with 

n 

P{y(t)-y(„)(t)>x}--5]e-"'^ 



n 
1=1 



In particular, for 1 < k < n 
k 



p{(y(t)-y(„)(t)) > (y(,)(t)-y(„)(t))} - i^exp(-a,(y(,,)(t)-y(„)(t))). 



i=l 

and we can use this equation to determine the shape of the capital distribu- 
tion curve. 

We wish to determine the shape of the graph of log k = log Yli=i e~"'^ , 
where log k is considered to be a function of x. For <p{x) := Y17=i e~"'^ , we 
have 

^ log = ^ log ^{x) = ^-^^ . 

For this derivative to be nonnegative, it suffices that the numerator be non- 
negative, and we have 

n n n 

i=l j = l i,j = l 



1 " 



2 
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SO the graph is convex. It follows that the capital distribution curve of the 



Example 8. To see different shapes of the expected capital distribution 
curve under different parameters apart from Examples [5] & [6l let us consider 
a pure hybrid market whose drift and volatility coefficients do not depend on 
ranks, except for the smallest (Atlas) stock. For example, take n = 5000 , 
gk = 0, 1 < k < n-1 , gn = c^i2n - 1) , 71 = -c* , ji = -2c^ , 2 <i <n, 
= 0.075, 1 < k < n , and pij = for 1 < i,j < n with a parameter 
c=K = 0.02. These parameters satisfy the assumptions of Proposition [3l We 
cannot apply Corollary[6l because the difference in (j5.2ip is positive on {p G 
S„ : p{k + 1) 7^ 1 } but negative on its (smaller) complement. The resulting 
expected capital distribution curve is convex; it is depicted in Figure [21 □ 

Example 9. Let us consider now a variant of this pure hybrid model, with 
a variance structure that is observed in practice. The parameters are the 
same as in Example [8l except for the different choices of the parameter c,, 
and for the rank-based variances ar := 0.075 + 6fcxl0"^ which are obtained 
from the smoothed annualized values for 1990-1999 data as in Section 5.4, 
page 109 of [12] (see page 2319 of [3j). The criterion from Corollary [6] cannot 
apply directly to this case, because the inequalities (I5.2ip do not hold for all 
p G S„ . The expected capital distribution curves under these parameters 
with (i) c=K = 0.02 , (ii) c^, = 0.03 , (iii) c^: = 0.04 are shown in Figure [3l 
The curve (i) is convex from the top rank to about the 25th rank, then turns 
concave until the lowest rank. The other curves (ii) & (iii) behave similarly. 



Example 10. Adopting the same parameter specifications in Example [9] 
(i) c* = 0.02 , except the rank-based drift, i.e., (iv) the upwind first ranked 
stock gi = -0.016 , gk = 0,2<k<n-l, g^ = (0.02)(2n - I) + 0.016 
and (v) the windward top 50 stocks gi = 52 = • • • = 550 = —0.016 , gt = 
0,51<fc<n — 1, Qn = (0.02)(2n — 1) -1-0.8 , we obtain concave curves as in 
Figure HI The observed average curve and the estimated curve of the first- 
order Atlas model for 1990-1999 (Figure 3 of [3j, page 2320) are concave. 
The statistical inference for the capital distribution curves is an interesting 
problem that we do not discuss here. □ 



Let us consider investing in the market of (j2.ip according to a portfolio 
rule n(-) = (ni(-), . . . , n„(-))' . This is an {.Ff}— adapted, locally square- 
integrable process with ^11=1 — 1 • Each Iii{t) represents the propor- 
tion of the portfolio's wealth V^{t) invested in stock i at time t, so 



pure hybrid market is likely to be convex. 



□ 



□ 



6. Portfolio Analysis 



(6.1) 
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For example, we may choose for every t G [0, oo) the vector of market 
weights Hi{t) , i = 1, . . . ,n as in (14. 2p . We shah call the resulting n(-) = 
//(•) the market portfolio. 

For a constant-proportion portfolio n(-) = vr G F" := {(tti, . . . , vr^)' G 
M" I Y17=i ^4 = 1} (which of course the market portfolio is not), the solution 
of (16. ip is given by 



(6.2) dlogVit) = J* {t)dt + ^TTid log Xi{t); 0<t<oo. 



i=l 



Here we denote by (aij(i))i<j •<„ = ^ 0^ i'^)) ^ 0^ i'^))' covariance process 
from p.4p . and introduce 

(6.3) 'y*{t) := -(^'^iTiaii{t) - 7riaij{t)7Tj^ ; < t < oo , 

i=l *ii=l 

the excess growth rate of the constant-proportion n(-) = vr G F". Thus, for 
a constant-proportion portfolio we can write the solution of ()6.ip . namely 
(6.4) 



y^(t) =u;.exp 



i=l '■^ ' i,j=l 



as in (2.4) of |T7], where Aij{-) = aij{t) dt ; we set A{-) := (^jj (•))i<i,j<n • 

6.1. Target Portfolio. Let us assume that, for every {t,uj) G [0,oo) x 
0, there exists a vector U*{t,u;) := (n^(t,cj), . . . ,n;(t,u;))' G F" that 
attains the maximum of the wealth V^{t,uj) over vectors tt G F" ; and 
that the resulting process !!*(•) defines a portfolio. Along with Cover [SJ & 
Jamshidian ^17j, we shall call this !!*(•) Target Portfolio, and 

(6.5) V^{t) := maxVit) , < t < oo 

the Target Performance for the model. (The quantity of (16. 5p is not nec- 
essarily equal to, and will typically be very different from, the performance 
V^*{-) of the portfolio U* .) 

The Target Performance Vi(-) exceeds the performance of the leading 
stock, of the value- line index (the geometric mean), and of any arithmetic 
average (such as the DJIA): to wit, taking Xi{0) = • • • = X„(0) = 1 , we 
have for every vector (qi, . . . , a„)' G F" := {(tti, . . . , 7r„)' G F"|7rj > , i = 
1, . . . , n} the almost sure comparisons 

n 1 / n " 

(6.6) K(-) > max [ max Xi(-) , ( J] Xj{-)) " , a,X,(- 

i=i i=i 

Under the assumptions of Proposition [H the limits Op of the average 
occupation times in (13. Sp exist almost surely, and so do the limits of the 
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average covariance rate := Huit^oo ^ij{T) /T , i.e., a°° := (a^) 
is 



dt 



lim 

T^oo 



1< j j<n 



pGS„ 



with Sp defined in (j2.3p . It follows from ()6.4p that the asymptotic long- 
term-average growth rate of a constant-proportion portfolio vr S F" is 

(6.7) ^hm 1 log F-(r) = 7 + ^ ( ^ vr,a,f - 7r,a,77r,) =: 7 + 7^ • 

i=l i,i=l 



Maximizing this expression over vr G F" amounts to maximizing over constant- 
proportion portfolios 7^ = (l/2)( X^iLi - TA,j=i ^i'^ij^j) ' the excess 
growth rate that corresponds to the asymptotic covariance structure. 

We shall call Asymptotic Target Portfolio a vector vf = (vfi, . . . , 7f„)' G 
that attains maxTrer" 7^ • We can regard this as an asymptotic growth- 
optimal portfolio over all constant-proportion portfolios, in the sense that 
liniT^oo (l/T) log{V''{T)/V*{T)) < holds a.s. for every vr G F" . 

Example 11. When there is no covariance structure by name, i.e., pij = 
for every 1 < i,j < n, we have Aij{-) = for i ^ j ■ In this case, we 
compute a target portfolio n*(-) as 



(6.^ 



mt) = (2A,,{t)Y^ 



-1 r 



1 ^nit) 



1 1 
+ 2+1^ 



n 



1 ^nit) 



log 



log 



X,{t) 
Xi{0) 



Xj{t) - 
,n. 



and an asymptotic target portfolio by 



(6.9) 



vr,; 



n 



lim U*it); 



, n, a.s. 



The portfolio vr has exactly the same long-term growth rate as the target 
performance in (|6.5p . namely 

1 /V^(T)\ 

(6.10) .1- r^°nT^)=0' 
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on the other hand, it outperforms the overall market rather significantly 
over long time horizons, namely 

lim ^ log ( ^^^^) = - y vfi(l - 7fi)a^ 

a.s., from the arithmetic mean - harmonic mean inequality. □ 

With Cover [8] and Jamshidian [17], we shall say that stock i is asymp- 
totically active, if for the expression of ()6.9p we have vf j > ; and that the 
entire market is asymptotically active, if all its stocks are asymptotically 
active, that is, if vf G F"^ := {(vri, . . . , 7r„)' G T" | vrj > , i = 1, . . . ,n} . 

Example 12. A sufficient condition for asymptotic activity of the model 
with n > 3 under the condition of Proposition [3l is obtained from (|6.9p as 

1 1 " 1 
(6.12) — — < -( — — ) , or equivalently 



n—l , ^ n n—1 

( E -^.(.) n \-])" < ;rb[E ( E n 

pGS„ i=l £=1 pGS„ j=l 

for every i = 1, • • • , n , with Apj defined in (j5.10p . This is the case in the 
constant variance model erf = • • • = cr^ . In general, it seems that the drift 
and volatility coefficients have non-trivial effects on the condition (j6.13p . □ 

6.2. Universal Portfolio. The Universal Portfolio Cover [8\ and Jamshid- 
ian [IT] is defined as 

J^„7rV^{t)dTT 

The wealth process of this portfolio is given by the "performance- weighting" 

V^t)d7T 

V^{t) = ^j — ; 0<t<oo, 

Jp„ dTT 

as can be checked easily. It follows from Theorem 2.4 of Jamshidian [T7] that 
the Universal Portfolio does not lag significantly behind the Target Portfolio: 
its performance lag is only polynomial in time under an asymptotically active 
model. To wit, there exists then a positive constant C , such that 

as T^oo 
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holds almost surely, thus also 

1 /V^(T)\ 1 /V^(T)\ 

lim -log (—4^) = lim -logf— ^) =0. 

In the context of the hybrid model, under the assumptions of Proposition 
[3] and of Example [T2l the Universal Portfolio attains the long-term growth 
rate of the asymptotic target portfolio vf (that is, the maximal long-term 
growth rate achievable through constant-proportion portfolios), as if the 
structure of the market were known. 

6.3. Growth-Optimal Portfolio. We shall call growth- optimal a portfolio 
tu(-) that satisfies the inequality Huit^oo {l/T) log{V^{T)/V'^{T)) < 
almost surely, for any portfolio n(-) . 

In order to find such a growth-optimal portfolio under no-name based 
correlation pij = for 1 < i,i < n , we need to maximize over vr G the 
quantity 



n ^ 1 



r2 

2 / J —1,1, \~ / I ) 

i=l " i=l 

where ji{t) = Epes„ l^^pC^W) 9p-i(j) +7^+7 is the i^^ element of 
G{Y{t)) of (j2.4p . By the Lagrange multiplier method, we obtain 

1 7i(t) +7(t) 

^iit) = 7T+ J ; ^ = l,...,n, 0<i<oo 

2 aii[t) 

where the constraint J2i=i = 1 is enforced by the multiplier 

n ^ -j^ ^ n ^^^^ 

The growth rate of this portfolio's performance V^{-) is 

1=1 1=1 1=1 

• In order to make some comparisons, let us specialize to the equal-variance 
case, i.e., cif = • • • = u,^ = ci^ with no name-based correlations pij = ; we 
obtain 

-I 2 1 1 " 

(6.14) ^,i„>^ - log y-(T) = 7 + :L (i _ _) + _ ( _ ^ 

k=l i=l 

On the other hand, the Universal Portfolio and the Asymptotic Target Port- 
folio have the same long-term growth rate, namely 

(6.15) hm i log V^T) = hm i log V^{T) = ^ + ^(l - - 
T— >oo 1 T— >oo 1 2 \ n 
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Under the conditions of ()2.2p and ()3.2p . we can verify 

n n 

(6.16) E^fe > 

k=l i=l 

To show (|6.16|) , we may assume without loss of generahty 71 > • • • > 7n 
and hence that there exists {61, ... , Sn-i)' S (M+)"~"'^ \ {0} such that gk = 

-{ik + 5k) for /c = 1, . . . , n - 1 , and gn = -7n + ((^i H h 5n-i) for ([2?2]) 

and (j3.2p . Then we obtain 

n n—1 

E^fc = E(^^ + '^^)' + (-^" + ('^1 + + 

fe=l i=l 

71 Tl — '1 71, — 1 77 

= E^' + E('^' + 25^(^^-7n)) + (E'^0' > 

i=l i=l i=l i=l 

Thus we observe from (j6.14p . (j6.15p . (j6.16p that the growth-optimal port- 
folio w{-) dominates in the long run both the universal portfolio n(-) and 
the asymptotic target portfolio vf , a.s. The advantage of the universal port- 
folio is that it can be constructed with total oblivion as to what the actual 
values of the parameters of the model might be; some of these may be quite 
hard to estimate in practice. By contrast, constructing the growth-optimal 
portfolio tu(-) requires knowledge of all the model parameters, and keeping 
track of the positions of all stocks in all ranks at all times. 



7. Appendix 

7.1. Proof of Lemma [Tl The stochastic exponential 
C(t) = exp [- [\^{u),dW{u))-\ f\Uu)fdu 



< t < 00 



is a continuous martingale, where ^{t) := 5 ^{Y{t))G{Y {t)) for < t < 00 
and := Yl^=i , x £ and {x, y) = Yl^=i ^jVj ' x,y G . Recall 

that S{-) , S~^{-) and G(-) in (|2.4p are bounded. By Girsanov's theorem 

W{t) := W{t) + / a~\Y{u))^x {Y{u))du, < t < 00 



is an n— dimensional Brownian motion under the new probability measure 
Q , locally equivalent to P , that satisfies 

(7.1) Q(C7) = E^(C(r) Ic) ; C G .Ft , < T < 00 . 
Thus, equation (j2.4p under P is reduced to 

(7.2) dY{t) = S{Yit))dW{t); 0<t<T, under Q . 



(7.3) A,*(t) = / l|,.(,)=o} -577V = 0; 0<t<cx,, l<5<2. 
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7.1.1. Local time of Bessel Process. Let us denote by .Z^(-) the square of the 
(5— dimensional Bessel process with < 5 < 2 . This process reaches the 
origin with probability one, but is subject to instantaneous reflection there, 
i.e., the local time ^z{') of Z{-) at the origin is zero: 

Kz{t) = 5 [ l{z(s)=a}ds =0 and so Leb({t : Z(t) = 0})=0 

holds for all t G [0, oo) almost surely. (Proposition XI. 1.5. of [21j) The 
same is true for the (5— dimensional Bessel process r''(-) = y^|Z(-)| with 
1 < 5 < 2, i.e., Leb({t : x^{t) = 0}) = almost surely, and hence 

5-1 ds 

V(^)=0}^) 

On the other hand, for 6 >2 the origin is never reached at all. We conclude 
that the local time of 5— dimensional Bessel process r''(-) is identically equal 
to zero, i.e., A.^s{-) = , for any 5 G (l,oo) . 

7.1.2. Comparisons with Bessel Processes. Now let us fix integers 1 < i < 
j < k < n . Under Q in (17. ip we shall compare the rank gap process 

r]{t) := max l£(t) — min Ym{t) 

£=i,j,k m=i,j,k 

with a Bessel process of dimension 6 > 1 , using Lemmata [5] and [6] below. 

We introduce the function g{y) := [{yi - y^Y + {yj - ykY + {yt - ViYY^"^ 
for y G , and note the comparison rj{-) > g(Y{-)) . An appli- 
cation of Ito's rule to g{Y{-)) yields the semimartingale decomposition 
dg{Y{t)) = h{Y{t))dt + dQ{t) , where we introduce the (n x 3) matrix 
Dijk '■= {di ,dj ,dk ) with (n x 1) vectors di := Cj — , dj := zj — , 
dk '■= ^k — , we denote by Cj , i = 1, . . . ,n the i*^ unit vector in , and 
{R{y) - l)Q{y) TT{D[^,S{y)S' {y)D,,,) 
''^'^■=^^) ' ''^'^■= ' 

y'DijkD.j^y 

(7.4) Z := {y G M" I g{y) = {y' D.^kD'^^f^y) = } , 







y=Y{s) 



(G)(t) 



/ Q{Y{s))ds] 0<t<oo. 
Jo 



Here note that under the assumption on ()2.3p . and because SDiji^D'-jj^ 
DijkD'ij^DijkD'-ji^ , we have 

y'DijkDl.^D.jkDl-^y pes„ £=i,...,n ^ 
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in M" \ Z ^ where A^^p,£ = are the eigenvalues of the positive- 

definite matrices SpSp for p G S„ , and so is strictly increasing when 

y(-) G \ Z . Now define the stopping time r„ := inf{t > | (e)(t) > u} , 
and note 

(S(n) := giJ{Tu)) = g{Ym + r ^(^W) + ^(^) ; < n < oo , 

JO 

where B{u) := @{tu) , < u < oo is a standard Brownian motion, by the 
Dambis-Dubins-Schwartz theorem of time-change for martingales. Thus, 
with d{u) := R(Y{tu)) we can write 

d&{u) = ^^"^^ ~ ^ du + dB{u) ; < u < oo , 0{O) = g{Y{0)) . 

The dynamics of the process are comparable to those of a Bessel 

process r'^(-) with dimension 5, generated by the same -B(-) and started 
at the same initial point g'(y(0)) . Since is positive definite un- 

der ([23]) and rank (Ajfc) = 2 , the (3 x 3) matrix L»^^5(-)5'(-)'Aifc is 
non-negative definite and the number of its non-zero eigenvalues is equal 
to rank(Z)^^5(-)5(-)'Ajfe) = 2. Let us denote by A^,p , i = 1,2,3 the 
eigenvalues of D'- j f^SpSp' Dijk for p G S„ . Then for R{-) in (|7.4p we obtain 

(7.6) R{-) > 6o := min [ ^^=i ^^lP | > i in M" \ ^ 

peE„ ymaxi<^<3 A^^p J 

and so ()(•) > 5o > 1 when Y(t.) G M" \ Z . By a comparison argument 
similar to that in the proof of Lemma 2.1 of [16], we may show that &{t) > 
x^o{t) for < t < oo a.s. Thus V3ri{t) > g{Y{t)) = (5((e)(i)) implies 
V3r]{t) > x^°{{@){t)) for < t < oo , a.s., and so we obtain the following 
result. 

Lemma 5. For the process Y{-) ofi^T§ with (KE), the multiple V3r]{-) of 
the rank-gap process dominates, a.s. under Q , a time-changed Bessel process 
r(-) '■= {{&) {■)) with dimension 6o as in (| 7. 6[) ; 



Q(\/37?(t) > r(t) ; < t < oo ) = 1 . 

Lemma 6. Under Q, the rank-gap process t]{-) satisfies {r]){t) < cit, 
< t < oo a.s. for some constant ci > and the local time Aj;(-) of ry(-) 
at the origin is identically equal to zero. 

In fact, since the diffusion coefficient matrix S{-) of Y{-) in (j7.2p is 
bounded and positive definite under (j2.3p . there exist such constant ci that 
{'n){t) < ci t for < t < oo a.s. Moreover, from (17. 5p and Lemma O there 
exists a constant C2 := minpgs„/=i,...,n ^i,p > , such that (6)(t) > C2 t for 
< t < oo a.s. It follows from the representation of local times (Theorem 
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VI. 1.7 of ^) and ([73]) with Lemma [5] that 

^r,{t) = — l{o<„(s)<.} d {r]){s) < lim ^ {o<Vzv{s)<e} d s 

<lim— / lrn<w.«w^\ as < hm / lrn<-.^^,^^^^ aw 



< \/3ciC2^ Ar^((0>(t)) = 0; 0<t<oo. 

Define an increasing family of events Ct '■= {Ay^{t) > for some t G 
[0,r]}, T > 0. By Lemma El we obtain Q(Coo) = and = Q(Q) = 
P(Q) for ^ > 1. Then P(A^(t) > for some i > 0) = P(U^iQ) = 
limi=^¥{C£) = 0. Thus the local time A^(t) of rank gap process ri{-) for 
{Yi{-) ,Yj{-) ,Yk{-) ) is zero for < t < oo a.s. under P. 

Since the choice of k is arbitrary, there is no local time generated by 
the rank gap process of any three coordinates. The rank gap process of more 
than three coordinates (e.g. maxi^h^ij^i^ y£(-) — minm=h,i,j,fc ^(O) dominates 
that of any three sub-coordinates. Therefore, by a similar argument as (|7.7p 
and its consequence, any local time of rank gap process of more than three 
coordinates is zero for < t < oo a.s. under P . 

To establish ()4.4p from this and (j4.3p . and thus complete the proof of 
Lemma [H consider any integers (ranks) l<a<i<m<b<n with 
b — a > 2 , and observe that we have almost surely: 



= A'^'^(t)= / l{z4s)=z,is)}d{Za{s)-Z,{s)) 
Jo 

^{Za(s)=Zt{s)} d{Za{s) - Zt{s)) + / l{Z,(,)=Zi,(.)} d{Zt{s) - Zm[s)) 



+ / l{Za{s)=Zi{s)} d{Zm{s) - Zb{s)) 

Jo 

Mz.is)=z,is)} d{A-'\s)+A''"^{s)+A^''is)) > f l{z4s)=z,(s)} dA^'™(^) 



/ 

and thus (fO]) as well. □ 

7.2. Proof of LemmaO For each A; = 1, . . . ,n — 1 the local time A'^''^"'"^(-) 
is a continuous additive functional of (H(-), *p.) with support in , and the 
expectation of A^''^'^'^{t) with respect to the invariant distribution •) is 
finite for t > . 



The a.s. equality Jg ^{Za{s)=z,,{s)} dA^''^{s) = follows readily from this, as 
does 

„j n k—1 

V.W>3}( E dA^^'{s)-Y,dA'^\s)) =0 
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It follows from the theory of additive functionals ^ that there is a finite 
measure on 5^ x S„ such that 

(7.8) ^eJ r g{E{s),^s)dA''''+Hs)] = \ I g{z,p) du^iz^p) 

for every bounded measurable function (7 : x S„ i-^ E . Let us denote by 
2^0fc(") = T^ki'i'^n) the marginal distribution on '^^ ■ The absolute continuity 
of ^Qk{') with respect to (n—1)— dimensional Lebesgue measure is argued by 
localization and the properties of Reflected Brownian motion as in Theorem 
7.1, Lemmata 7.7 and 7.9 of [Hj. 

Now, by an application of Ito's rule, for / G C^((]R+)"^^) we obtain 

/(H(r)) = /(H(0)) + f {Vf{E{s)),dr^^\s)) 
Jo 

+ [T^kf]{^{s))dk^'^+\s)+ f [Af]{E{s),^,)ds; T>0, 

where is the martingale part of C,{ ) and and A are differential 

operators defined in (15. 3p . Taking expectations with respect to P and then 
integrating for the initial values with respect to the stationary distribution 
•) with Fubini's theorem and (j7.8p . we obtain 

= 9 E / i'^'^f] ^ ^Ofc {z)+T [Af] {z, p)d u{z, p) . 



k=l 



Dividing by T > , we obtain the basic adjoint relationship ()5.6p . 

7.3. A Sanity Check of Corollary [4l In this section we verify that 
{dk,i)i<i,k<n in (|5.16p Satisfy (|4.10p . Since O/^^i is homogeneous in the prod- 
uct 11^=1 [-4(fTj+cT|_^;^)"^] , it suffices to show X]fc=i Gk,i{gk+li) = where 
we use the modifications O^^i := J2{p(k)=i} ' 

n— 1 n—1 j 

- ( E n ~^^^^) ' n ' ^p.^- ^= E + ^pw) 

qGS„ j = l j = l 1=1 

of {0k^i,6p, Apj) , 1 < i, j, k < n , p G S„ for notational simplicity. Note 
that Xp^n = from (j2.2p for p G S„ . 

First, observe for I = 2, . . . ,n and i = 1, . . . , n , 

(7.9) Yl v-i^P+ E (9i+7iK= E v^p- 

{p:p(£-l)=i} {p:pW=i} {P-P{i)=i} 

In fact, for every i,i define another permutation p from a (fixed) permu- 
tation p G {q G S„ : q{£ — 1) = i} by 

p(A;), A; = l,...,^-2, ^ + l,...,n, 
pik) ■.= pik;p) = { p{i), k=l-l, 
% , k — ^ . 
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which is obtained by exchanging (i — l)st and £th elements of p G {q S 
S„ : q(^— 1) = i} , and also define M := (X^qg^n 11^=1 -^qj) ""^ here. Then 
Xpj = Xpj for j ^ £ — 1 and hence the left-hand of ()7.9p is 

n—1 n—1 
{p:p(^-l)=j} j = l {p:p{i)=i} i = l 



n—1 n—1 

a; . 
'J 



E ^ n ^pi+ E i9i+imwU~^p] 

{p:p{l)=i} {p--W)=i} i=l 

n-1 

E [^p/-i + 5^ + ^p(^)] • ^ n ^p,i = E ^P'^^p ' 

{p:pW=*} J = l {p:pW=i} 



which is the right-hand of (j7.9p . Now applying (j7.9p for £ = 2, . . . , n , we 
obtain 

n n 

E(gfe + 7i)^fc,i = (91 + 7i)6'i,j + (92 + 7j)^2,j + E'^^'^ + li)^k,i 

k=l k=3 

n 

= E -^p.2^p + E^^'^ + iiWk,i = ■■■ = E Ap,„6ip = , 

{p:p(2)=i} fc=3 {p:p(n)=j} 

for i = 1, . . . , n , because Ap_„ = for p G S„ . Therefore, (j4.10p is satisfied. 
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Figure 1. Different values of 
{Qk,i\ for (/c, i) , when the pa- 
rameters are specified for an 
extreme case in Example HI 
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Figure 2. Expected capital 
distribution curve for the pure 
hybrid model in Example [HI 




Figure 3. Expected capital 
distribution curves for the hy- 
brid model in Example [9j 



Figure 4. Expected capital 
distribution curves for the hy- 
brid model in Example llOi 



